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Resumen 


Los flujos verticales, dentro del ciclo hidrológico, son una de las variables 
de mayor relevancia en zona de llanura, dado que las pendientes varían 
entre O y 5%, y los flujos horizontales no son significativos. En este 
sentido, la evapotranspiración juega un rol fundamental en el manejo 
hídrico, ya que alrededor del 85% del agua que sale del sistema lo hace 
mediante este proceso, requiriendo una cuantificación precisa. El objetivo 
principal de este trabajo es calcular la evapotranspiración potencial y real 
(ETp y ETr) con datos de satélite y reanálisis mediante el uso de la 
plataforma Google Earth Engine. Para su cálculo se tomó la ecuación de 
Priestley-Taylor (PT) combinado con información de humedad de suelo, 
con una resolución espacio-temporal de 250 x 250 m cada ocho días, en 
la región pampeana argentina (RPA). El producto se valoró en siete 
estaciones de la RPA, cuyos resultados mostraron, para la ETp, un R- 
RMSE (error cuadrático medio robusto) de 0.5 mm d*; un error 
sistemático (mediana) de 0.3 mm di, y el error aleatorio (RSD- 
desviación estándar robusta) de 0.5 mm d**; mientras que para la ETr, 
estos valores son de 0.6, -0.2 y 0.5 mm d*?!, respectivamente. Los 
resultados globales muestran que el método utilizado es una herramienta 
válida para caracterizar la ET en la RPA y que se puede usar para analizar 
su variabilidad espacio-temporal en diferentes condiciones extremas y 


realizar estudios ambientales aplicados. 


Palabras clave: MODIS, NOAA, SMAP, Priestley-Taylor, región 


pampeana argentina. 
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Abstract 


Vertical flows, within the hydrological cycle, are one of the most relevant 
variables in the plains, since slopes vary between O and 5%, and 
horizontal flows are not significant. In this sense, evapotranspiration plays 
a fundamental role in water management since about 85% of the water 
leaving the system does so through this process, requiring precise 
quantification. The main objective is to calculate potential and actual 
evapotranspiration (ETp and ETa) with satellite and reanalysis data using 
Google Earth Engine platform. For its calculation, the Priestley-Taylor (PT) 
equation combined with soil moisture information was used, with a spatio- 
temporal resolution of 250 x 250 m every eight days, in the Argentine 
Pampas region (APR). The product was valued in seven stations of the 
APR, whose results showed, for ETp, an R-RMSE (Robust Root Mean 
Square Error) of 0.5 mm d*, a systematic error (Median) of 0.3 mm d”*!, 
and the random error (RSD- Robust Standard Deviation) of 0.5 mm d”!; 
while, for ETa, these values are 0.6, -0.2 and 0.5 mm d”!, respectively. 
The overall results show that the method used is a valid tool to 
characterize ET in the APR and that it can be used to analyze its spatio- 
temporal variability under different extreme conditions and to carry out 


applied environmental studies. 


Keywords: MODIS, NOAA, SMAP, Priestley-Taylor, Argentina Pampas 


region. 
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Introducción 


Los flujos verticales dentro del ciclo hidrológico tienen un rol fundamental, 
en particular la evapotranspiración (ET), ya que en zonas de llanuras se 
pierde alrededor del 85% del agua del sistema mediante este proceso 
(Weinzettel 8 Usunoff, 2001). 


En este sentido, y dado que dentro de los objetivos de la agenda 
2030 para el Desarrollo Sostenible de las Naciones Unidas se destaca la 
importancia de la disponibilidad y gestión sostenible del agua y 
saneamiento, así como promover la agricultura sostenible y tomar 
medidas urgentes en torno a combatir el cambio climático y sus impactos, 
el análisis y estudio del ciclo hidrológico, y en particular de los flujos 
verticales, en los cuales la atmósfera interactúa con el agua superficial, 
las superficies con vegetación y los suelos desnudos, es de vital 
importancia para el cuidado y gestión de los recursos tanto a nivel local 
como global (McMahon, Peel, Lowe, Srikanthan, € McVicar, 2013; Liu, 
Yao, £ Wang, 2019; Peng et al., 2019; Tasumi, 2019; Xiang, Li, Horton, 
8: Feng, 2020; Degano, 2021a). 


Con respecto a la evapotranspiración (ET), cuantificarla espacial y 
temporalmente tiene importancia en la gestión de los recursos hídricos; 
en la determinación de la huella hídrica (Liu € Yang, 2010; Olivera- 
Rodríguez et a/., 2020), y en la mejora de la productividad del agua y la 
seguridad alimentaria (Garg et a/., 2012; Liu, Zehnder, € Yang, 2009). 


Según la escala de trabajo que se requiera, la ET se puede estimar 


con diversos métodos: (1) a escala puntual: porómetros, medidores de 
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flujo de savia en la planta y lisímetros de pesada (medida directa) (Allen, 
Pereira, Howell, 8 Jensen, 2011; Silicani, 2015; Faramiñán, Carmona, 
Rivas, Degano, 8l Olivera-Rodríguez, 2021a; Pendiuk, 2022); (2) a escala 
de parcela: balances de agua, relación de Bowen, sistemas Eddy 
Covariance y centillómetros, entre otros; (3) a escala regional: datos de 
satélite (Rivas 8 Carmona, 2010; Ocampo € Rivas, 2013; Degano, 2021a; 
Degano et al., 2021b; Degano, Rivas, Carmona, Niclós, € Sánchez, 
2021c). 


Estos métodos presentan errores asociados que pueden estar 
relacionados con el sensor (sesgo de calibración, funcionamiento, 
operación co ubicación incorrecta); con la interpretación o el 
procesamiento erróneo de los datos; con las características de la 
vegetación; con la medición, que incluye componentes aleatorios 
asociados con la resolución de las lecturas del sensor, el ruido electrónico, 
el ruido inducido mecánicamente, las respuestas térmicas de los sensores, 
la gestión de la vegetación y el agua del suelo, así como con otros errores 
aleatorios específicos para el tipo de sistema de medición (Allen et al., 
2011; Gao et al., 2020). 


Es importante destacar también los errores relacionados con los 
diferentes métodos de cálculo de ET (Allen et a/., 2011), en particular el 
balance de agua en el suelo, con un error típico de 10 a 30%; el balance 
de energía con teledetección, cuyo error oscila entre 10 y 20%; los 
gravímetros superconductores, que disminuyen el error del método a 
valores de 3 a 4% (Pendiuk, 2022). 


En general, la medición directa de la ET ofrece grandes dificultades 


tanto a escala regional como global, dado que los factores que intervienen 
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en este proceso son diversos, y variables en el tiempo y en el espacio; es 
donde adquieren mayor relevancia los datos de satélite, que facilitan la 
determinación de la variación espacial y temporal de los diferentes 
términos de la ecuación de balance de energía (Brutsaert, 1984) a escalas 
espaciales que pueden ir desde parcelas hasta escala de cuenca, así como 
la obtención de datos a nivel global, al facilitar el estudio y análisis de la 


variabilidad climática. 


Dentro de las plataformas de mayor envergadura e importancia en 
cuanto a las variables biofísicas se encuentran distintos satélites que 
portan sensores a diferentes resoluciones espectrales, espaciales y 


temporales (Tabla 1). 
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Tabla 1. Características de las distintas misiones de satélite. 


nd y : . Resolución Resolución 
Misión/sensor Años activo Procedencia A € 
espacial (m) | temporal (días) 
Landsat 1-9 (Tatem, Goetz, 8 Hay, 2008) | 1972-Actualidad | NASA!/USGS? 30-60-100 16 
CERES? 
1997-Actualidad NASA/JAXA* 10 1 


(Smith et a/., 2011) 


MODIS? (Terra y Aqua) 


Cutie ebal 1998) 1999-Actualidad NASA 250-500-1 000 pl 
Copernicus Sentinel 1-6 2014-Actualidad ESA? 10-20-60 10 

SMOS” (Barré, Duesmanmn 8, Kerr, 2008) | 2009-Actualidad ESA 50 000 3 
SMAP? (Entekhabi et al., 2010) 2015-Actualidad NASA 36000 2-3 
SAOCOW? 1 A y B (Giraldez, 2004) 2018-Actualidad | CONAE*%YASI!! 10-100 16 


National Aeronautics and Space Administration. 
2United States Geological Survey. 

3Clouds and the Earth's Radiant Energy System. 
“Japan Aerospace Exploration Agency. 

Moderate Resolution Imaging Spectroradiometer. 
SEuropean Space Agency. 

7Soil Moisture and Ocean Salinity. 

8Soil Moisture Active Passive. 

“Satélite Argentino de Observación Con Microondas. 


10Comisión Nacional de Actividades Espaciales. 


liAgencia Espacial Italiana. 
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Como complemento a la información que brindan las misiones 


satelitales, es importante destacar los conjuntos de datos de reanálisis 


(Tabla 2). Se trata de bases de datos que permiten desarrollar un registro 


completo de la variabilidad climática en el tiempo. Tienen en cuenta las 


observaciones y un modelo numérico que simula uno o más aspectos del 


sistema terrestre, y se combinan objetivamente para generar una 


estimación sintetizada del estado del sistema. Por lo general se extienden 


a lo largo de varias décadas y cubren todo el globo (Wielicki et a/., 1996; 
Rienecker et a/., 2011; Poli et al., 2016). 


Tabla 2. Características de los productos de reanálisis. 


a . Resolución | Resolución 
Base de datos de reanálisis Periodo Procedencia : 

espacial temporal 

ERA5 (Muñoz-Sabater, 2019) 1950-Actualidad ECMWF! 3 000 m 3 horas 

JRA-55? (Kobayashi 8: National Center for 
1958-Actualidad JMA? 1,250 6 horas 
Atmospheric Research Staff, 2015) 
NOAA*-CFSV2* 
1979-Actualidad NCEP?/CFS 22 264 m 6 horas 
(Saha et al., 2011) 
MERRA 2-GMAO” (Gelaro et al., 2017) 1980-Actualidad NASA 50 km Diario 


“European Centre for Medium-Range Weather Forecasts. 


2Japanese 55-year Reanalysis. 


3Japan Meteorological Agency. 


“National Oceanic and Atmospheric Administration. 


Climate Forecast System versión 2. 


ENational Centers for Environmental Prediction. 


7Modern Era Retrospective-analysis of Research and Applications of Global Modelling 


and Assimilation Office. 
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Las ventajas que presentan las plataformas satelitales y las bases 
de datos de reanálisis dieron lugar a un gran número de aplicaciones que 
permiten obtener variables de superficie a partir de datos registrados en 
diferentes sectores del espectro electromagnético. En relación con esto, 
es de gran utilidad la plataforma Google Earth Engine (GEE), que combina 
diferentes catálogos de imágenes satelitales y conjuntos de datos 


geoespaciales en un mismo sitio, y opera integramente en la nube. 


La plataforma GEE de uso gratuito proporciona acceso a (1) 
imágenes de teledetección disponibles públicamente y productos 
procesados para su uso directo con una aplicación web de exploración; 
(2) procesamiento paralelo de alta velocidad y algoritmos de aprendizaje 
automático que utilizan la infraestructura informática de Google, y (3) una 
biblioteca de interfaces de programación de aplicaciones (API) con 
entornos de desarrollo que admiten lenguajes como JavaScript y Python. 
Juntas, estas características principales permiten a los usuarios descubrir, 
analizar y visualizar grandes datos geoespaciales sin necesidad de acceso 
a súpercomputadoras o experiencia en codificación especializada. El 
desarrollo de GEE ha generado entusiasmo y compromiso en los campos 
de la ciencia de datos geoespaciales y la teledetección. Sin embargo, 
después de una década desde que se lanzó GEE, su impacto en la ciencia 
geoespacial y la teledetección no se ha explorado a profundidad (Mutanga 
8: Kumar, 2019; Tamiminia et a/., 2020). 


En este sentido, el objetivo general del presente trabajo es estimar 
la evapotranspiración potencial (ETp) y real (ETr) aplicando la ecuación 


de Priestley-Taylor combinada con datos de humedad de suelo (Hs), 
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utilizando la plataforma GEE, con una resolución espacial de 250 x 250 m 
y temporal de ocho días. Para este fin se incorporan datos de entrada a 
la ecuación de ET de NOAA/CFSV2, MODIS y SMAP, que serán evaluados 
con la base de datos NASA Power y con datos de campo. Además, se 
valorará el producto de ETp y ETr con datos de ET medidos in situ. Por 
último, se busca analizar el comportamiento de la ET a escala espacial por 


medio de mapas regionales. 


Metodología 


Desarrollo del modelo. Ecuación de Priestley-Taylor 


Haciendo uso de la plataforma en línea GEE se determinó la ETp y la ETr 
con una resolución espacial de 250 x 250 m y temporal de ocho días, este 
rango temporal definido es capaz de representar el desarrollo vegetativo 
de las plantas/cultivos. Para el cálculo de ETp (en este caso es sinónimo 
de ETo, porque se trabaja sobre parcelas experimentales de pasto corto), 
se utilizó como base la ecuación de PT (Ecuación (1)) (Priestley € Taylor, 
1972) para obtener un modelo robusto, sencillo y de fácil aplicación. La 
consistencia del método está demostrada por diferentes autores 
(Stannard, 1993; Fisher, Tu, € Baldocchi, 2008; Rivas € Carmona, 2010; 
Miralles et a/., 2011; Marshall, Tu, 8 Andreo, 2020; Chandrasekar et al., 


2022; entre otros), que utilizan la ecuación de PT para determinar ET: 
ETppr = a(55) (Rn — G) (1) 
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siendo: 
a = equivalente a 1.26 
A = pendiente de la curva de saturación de vapor (Ecuación (2)) (kPa 0C- 
y 
y = constante psicrométrica (Ecuación 4) 
Rn = radiación neta (Ecuación 6) 


G = flujo de calor en el suelo (este parámetro se desprecia en proporción 


al valor de Rn (Rivas y Ocampo, 2009)). 


_ 4098€s (2) 


— (Ta+237.3)? 


donde: 


es = equivale al promedio de la presión de saturación de vapor del aire 
(kPa) 


Ta = temperatura del aire (0C) 


e, = 0.6108 exp (==) (3) 


Ta+237.3 


siendo: 


exp = 2.7183 (base de In) elevado a una potencia 
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pl CpPatm 
WE (4) 


donde: 
Cp = calor específico del aire y equivale a 1.013x10"3 (MJ kg"? 0C*1) 
Patm = presión atmosférica (kPa) 
£ =cociente del peso molecular del vapor de agua/aire seco = 0.622 


Á = equivale al calor latente de vaporización (MJ kg"*) (Ecuación (5)) 


2 = 2.501 — (0.002631 Ta) (5) 

Rn = (1 alb)Rs* + (e RI! — eo Ts?) (6) 
donde: 

alb = albedo (adimensional) 


Rs* = radiación de onda corta entrante (Wm-2) 


e = emisividad de la superficie (adimensional) (Ecuación (7); Van der 
Griend € Owe, 1993) 


RI* = radiación de onda larga entrante (Wm-2) 
v = constante de Stefan-Boltzmann 5.67 x 107 (Wm-? K-*) 
Ts = temperatura superficial (K) 


(eoTs*) = equivale a la RI? (radiación de onda larga saliente de la 


superficie) 
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€ = 1.0094 + 0.047 x In NDVI (7) 


donde: 
NDVI = índice de vegetación de diferencia normalizada (adimensional) 


Por otro lado, teniendo en cuenta el método de PT, la ETr 
corresponde a la relación entre la ETp estimada con PT y la fracción 
evaporativa (FE) (Ecuación (8); Davies € Allen, 1973): 


ETr =FE ETppr (8) 
donde: 
FE (Ecuación (9)) = modificador de las condiciones de humedad, 


relacionando la demanda máxima de agua de la atmósfera con la 
disponibilidad de agua real del sistema suelo-planta (Rivas € Ocampo, 
2009) 


FE = HS¡— HSmin (9) 


HSmax — HSmin 


donde: 


Hs = humedad de suelo 


Hs; = Hs instantánea 
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HSmin = Hs mínima o punto de marchitez 
HSmax = Hs máxima o capacidad de campo 


La estimación de este parámetro permite conocer el agua disponible 
en el suelo, que da lugar al proceso de ETr. Varía en particular con las 
características del suelo, su textura y estructura, así como por el tipo de 


suelo. 


Datos de satélite 


Dado que en la plataforma GEE se encuentra una gran cantidad de 
catálogos de datos de satélite, al igual que reanálisis de manera gratuita 
y disponible para las unidades académicas, se calcularon la ETp y ETr 
utilizando las bases de datos disponibles en esta plataforma (Tabla 3), 


según las ecuaciones descritas en el apartado anterior. 
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Tabla 3. Catálogos de GEE utilizados como datos de entrada para la 


ecuación de PT. 


E Resolución 
; Resolución 
Variable Catálogo espacial Unidad Periodo 
temporal 
(m) 
Ta: temperatura del 1979-01-01 
NOAA/CFSV2/FOR6H 6 horas 22 264 K 
aire Actualidad 
Patm: presión 1979-01-01 
a NOAA/CFSV2/FOR6H 6 horas 22 264 Pa 
atmosférica Actualidad 
Rs!: radiación de onda 1979-01-01 
NOAA/CFSV2/FOR6H 6 horas 22 264 Wm-? 
corta entrante Actualidad 
RI!: radiación de onda 1979-01-01 
NOAA/CFSV2/FOR6H 6 horas 22 264 Wm-? 
larga entrante Actualidad 
Ts: temperatura de ; 2000-02-18 
MODIS/061/MOD11A2 8 días 1 000 K 
superficie Actualidad 
, 2000-02-18 
NDVI MODIS/061/MODO9Q1 8 días 250 Adim. 
Actualidad 
2000-02-24 
Albedo MODIS/006/MCD43A3 Diario 500 Adim. 
Actualidad 
NASA_USDA/HSL/ 
, 2015-04-02 
Humedad de suelo SMAP10KM_ 3 días 10 000 mm 
Actualidad 


soil_moisture 


El catálogo NOAA/CFSV2/FOR6H disponible en GEE proviene del 


National Centers for Environmental Prediction (NCEP) Climate Forecast 


System (CFS), que es un modelo que representa la interacción entre la 


atmósfera, los océanos, la tierra y el hielo marino. CFS fue desarrollado 
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en el Environmental Modeling Center (EMC) en NCEP, donde se creó el 
Climate Forecast System Reanalysis (CFSR), que desarrollaron el 
conjunto de datos CFSv2 (Kalnay et al., 1996; Saha et al., 2011). 
Presentan datos de variables hidrometeorológicas promedio cada seis 
horas (0000, 0600, 1200 y 1800 UTC). En particular, la variable Ta que 
se encuentra en este catálogo es la temperatura a una altura de 2 m. 
Respecto de la Patm, la Rs! y la Rl', se obtienen valores promedio en seis 


horas. 


La variable Ts se extrajo a partir del catálogo 
MODIS/061/MOD11A2. El producto MOD11A2 versión 6 proporciona una 
temperatura de la superficie terrestre promedio de ocho días en una 
cuadrícula de 1 200 x 1 200 kilómetros. Cada valor de píxel en MOD11A2 
es un promedio simple de todos los píxeles Ts MOD11A1 correspondientes 
recopilados dentro de ese periodo de ocho días (Wan, 2019). Es 
importante indicar que con este producto se obtiene el valor de la Ts 
diurna y nocturna; para la aplicación de este producto se realizó un 


promedio de la Ts para obtener un valor diario. 


Los valores de NDVI se obtienen haciendo el cálculo según la 
Ecuación (10), empleando como entrada el catálogo de datos 
MODIS/061/MODO9Q1, que proporciona datos de reflectancia espectral 
de la superficie de las bandas 1 y 2, con una resolución de 250 x 250 m, 
corregida para condiciones atmosféricas (gases, aerosoles y dispersión de 
Rayleigh). Junto con las dos bandas de reflectancia, también se incluye 
una capa de calidad. Para cada píxel se selecciona un valor de todas las 
adquisiciones dentro del compuesto de ocho días en función de la alta 


cobertura de observación, el ángulo de visión bajo, la ausencia de nubes 
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o sombra de nubes, y la carga de aerosoles (Vermote, Roger, 8: Ray, 
2015). 


_ (NIR=Red) 
AS (NIR+Red) (10) 


siendo: 
NIR = banda de reflectividad del infrarrojo cercano (sur_refl_b02 en GEE) 
Red = banda perteneciente reflectividad en el rojo (sur_refl_b01 en GEE). 


El conjunto de datos del modelo MCD43A3 versión 6 de albedo es 
un producto diario de 16 días. Proporciona reflectancia hemisférica 
direccional (albedo de cielo negro) y reflectancia bihemisférica (albedo de 
cielo blanco) para cada una de las bandas de reflectancia de superficie 
MODIS (1 a 7), así como para las bandas del visible, infrarrojo cercano y 
onda corta. Cada imagen diaria de 500 x 500 m se genera utilizando 16 
días de datos, centrados en el día de interés. También se proporciona una 
banda de calidad para cada una de las 10 bandas de albedo (Strahler et 
al., 1999). 


Los datos de Hs global SMAP mejorado de NASA-USDA proporcionan 
información sobre la Hs en todo el mundo, con una resolución espacial de 
10 km. Este conjunto de datos incluye Hs superficial y subterráneo (mm), 


perfil de Hs (%), y anomalías de Hs superficial y subterráneo (-). 


El conjunto de datos se genera mediante la integración de 
observaciones de Hs de nivel 3 de SMAP, las anomalías de Hs se 


calcularon a partir de la climatología del día de interés. La climatología se 
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estimó con base en el registro completo de datos de la observación 
satelital SMAP y el enfoque de ventana móvil centrada en 31 días. La 
asimilación de las observaciones de Hs SMAP ayuda a mejorar las 
predicciones basadas en modelos, en particular en zonas que carecen de 
datos de precipitación de buena calidad (O'Neill, Chan, Njoku, Jackson, € 
Bindlish, 2016). Este conjunto de datos fue desarrollado por el 
Hydrological Science Laboratory at NASA's Goddard Space Flight Center 
en cooperación con USDA Foreign Agricultural Services y USDA Hydrology 


and Remote Sensing Lab. 


Procesamiento y valoración de los datos de satélite 


La secuencia de procesamiento aplicada en lenguaje JavaScript de la 
plataforma GEE integra la definición del sitio de interés; aplicación de 
funciones de filtrado espacial y temporal de los catálogos para la zona de 
estudio (RPA); selección y conversión de datos de variables involucradas 
en el cómputo, y desarrollo del algoritmo de PT para los valores de ETp 
con disgregación del dato a 250 m de resolución espacial. Este resultado 
se multiplica por la FE (que tiene como datos de entrada a la Hs), para 
obtener la ETr a la resolución deseada. Dicha resolución se logra aplicando 
el remuestreo con interpolación bilineal para cada conjunto de datos 
(Smith, 1981). 


Este método de remuestreo determina el nuevo valor de píxel 
basándose en el promedio de distancia ponderada de los cuatro centros 
de píxeles más cercanos. Es válido para zonas de llanura, donde el 


procesamiento del relieve es omitido y otorga resultados esperables para 
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mejorar la resolución espacial del producto (Salas-Aguilar € Paz-Pellat, 
2018) (Figura 1). 


» Definición del sitio de interés 

» Declaración de la función de 
corte (RPA) 

» Declaración de la paleta de 
colores 

* Fitrado temporal 


IT » Obtención de los parámetros 
Albedo NDVI 
(MODIS/006/MCD43A3) (MODIS/061/MODO09Q1) 
* Aplicación de la Ecuación (1) 
* Remuestreo a 250 m 


ETp (250 m) 


Ta, Patm, Rs, RI Ts 
NOAA/CFSV2/FOR6H (MODIS/061/MOD11A2) 


FE (Hs) 
(NASA_USDA/HSL/ 
SMAP10KM_ 
soil_moisture) 


I + Aplicación de la Ecuación (8) 
ETr (250 m) 


Figura 1. Diagrama de flujo del procesamiento de los datos de satélite 


para el cálculo de la ETp y ETr. 
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Para la valoración de los datos, en primera instancia se evaluaron 
particularmente las variables de entrada para estimar ET (Rs!, RlI*, Ta y 
Ts), que son los parámetros que aportan en el balance de energía, y por 
consecuencia, en el balance hídrico (principalmente la Rs'; Faramiñán et 
al., 2021a). Para llevar a cabo tal evaluación, se compararon las variables 
de satélite obtenidas de GEE (NOAA/CFSV2 y MODIS), con los obtenidos 
de NASA Power (base de datos validada en la RPA por Orte et al., 2021; 
Faramiñán, Degano, Carmona, €: Rodríguez, 2021b; Garcia-Bu-Bucogen, 
Piccolo, 8 Bohn, 2022), cuyos datos provienen de reanálisis (CERES, 
MERRA2), en el periodo 2015-2021. 


Por un lado, Orte et al. (2021) compararon los datos diarios de Rs 
de NASA-Power con mediciones en tierra en las ocho estaciones de la red 
de irradiancia solar Saver-Net instalada en el sur de América del Sur. 
Realizaron un análisis de regresión lineal para analizar la concordancia 
entre los datos satelitales y las mediciones terrestres, donde el coeficiente 
de determinación mostró una correspondencia relevante, con un valor 
medio de 0.95. El error absoluto medio (MAE) (0.63 kWh m-? d1) y el 
RMSE (0.48 kWh m2 d:1) reflejan una diferencia esperable. 


Por su parte, Faramiñán et al. (2021b) obtuvieron el valor de ETr 
utilizando NASA-Power como datos de entrada del algoritmo Support 
Vector Machine (SVR) en la RPA. Compararon los resultados obtenidos 
con valores de ETr (balance hídrico) con base en información de 12 
estaciones agrometeorológicas de la ORA (1983-2012). Observaron 
errores estadísticos medios de 0.39 + 0.07 mm d*!, 0.54 + 0.09 mm d*! 


y 0.67 + 0.07 para MAE, RMSE y R?, respectivamente. Los resultados 
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muestran la factibilidad de aplicar algoritmos de aprendizaje automático 


para la obtención de valores de ETr en planicies agrícolas. 


En este sentido, se considera que el conjunto de datos NASA Power 
es consistente para analizar los datos de satélite (GEE) en la RPA. 
Asimismo, con el fin de constatar la información con datos de campo, se 
verificó la respuesta de las variables de entrada con datos medidos in situ, 
en una parcela de cebada ubicada en las cercanías de la localidad de 
Tandil (-37.49; -58.90) durante el periodo julio-diciembre de 2019 


(sensores y variables detalladas en Faramiñán et al., 2022). 


Por otro lado, para evaluar el producto final de ETp y ETr, se 
obtuvieron los datos de entrada con una resolución temporal promedio de 
ocho días, mientras que para la realización de los mapas se optó por una 
resolución temporal anual para analizar los cambios en ese periodo de 
tiempo. La valoración del producto de ETp y ETr se hizo en el periodo 
2015-2021 debido a la disponibilidad de datos (tanto de campo como de 
satélite), y a la variabilidad en las precipitaciones, abarcando temporadas 


de máximos, mínimos y promedio anuales. 


Datos in situ, zona de estudio 


Los datos locales los determinó la Oficina de Riesgo Agropecuario (ORA). 


Calculan la ETo con datos meteorológicos medidos en estaciones por el 


Servicio Meteorológico Nacional de Argentina (SMN), y se basan en el 
método de FAO-Penman Monteith (PM) (Allen, Pereira, Raes, 8 Smith, 
1998; Monteith € Unsworth, 1990). 
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Por otro lado, los valores locales de ETr se obtienen realizando un 
balance hidrológico, de resolución temporal diaria, donde se tienen en 
cuenta los ingresos y egresos de agua del sistema, despejando el término 
de ETr para su cálculo (Allen, et a/., 1998; Basualdo, 2011). 


En cuanto a la RPA, es una región de planicie de alrededor de 520 
000 km? de superficie, que se encuentra en el centro-este de la República 
Argentina. En su extensión, incluye en parte cinco provincias agrícolas: 
sur de Entre Ríos, sureste de Córdoba, sur de Santa Fe, este de La Pampa 
y la provincia de Buenos Aires (exceptuando el extremo sur). Una estación 
meteorológica de referencia en zonas de llanuras representa un área de 
cobertura de 127 km de radio, delimitada por los círculos rojos que se 


encuentran en la Figura 2 (según la WMO, 1994; Degano et a/., 2021c). 
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Figura 2. Ubicación de la RPA y de las estaciones de validación (1-7). 
Distribución espacial del área de influencia de las estaciones 
hidrometeorológicas (círculo rojo) e isolíneas de ETr (INA, 2002) 


(modificado de Degano et al., 2021c). 
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Asimismo, la RPA se subdivide en 12 ecorregiones (Pereyra, 2003): 
planicies poligenéticas de entre ríos (D25); pampa pedemontana 
cordobesa (D26); pampa arenosa (D27); pampa endorreica (D28); delta 
del Paraná (D29); pampa ondulada (D30); pampa deprimida (D31); 
planicies litorales pampeanas (D32); sierras septentrionales bonaerenses 
(D33); pampa interserrana (D34); sierras australes bonaerenses (D35), 
y planicies poligenéticas de la pampa (D36), donde el clima general es 
templado-cálido y en particular varía de húmedo/subhúmedo (estaciones 
1. Concordia; 2. Rosario; 3. Laboulaye; 4. Junín, y 6. Tandil, que 
constituyen la pampa subhúmeda-húmeda) a semiárido (estaciones 5. 
General Pico y 7. Pigúé, que determinan la pampa semiárida) (Kottek, 
Grieser, Beck, Rudolf, 8 Rubel, 2006). 


En la Figura 2 se muestran además isolíneas de ETr definidas por el 
Instituto Nacional del Agua (INA, 2002). Ello determina la variabilidad de 
este parámetro, dependiendo de la zona de la RPA en la que se hace foco. 
Es mayor en las zonas más húmedas (noreste de la RPA), y disminuye 
hacia las zonas semiáridas de la región, variando entre 1 050 y 400 mm 


por año. 


Análisis estadístico 


Para analizar el comportamiento del producto de ETp y ETr obtenido es 
necesario evaluar su respuesta, para ello se utilizaron distintas métricas 
estadísticas Me (mediana), para analizar el error sistemático; RSD 


(desviación estándar robusta (Ecuación (11); (Wilrich, 2007), para 
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corroborar el error aleatorio; R-RMSE (error cuadrático medio robusto - 
Ecuación (12)); R? (coeficiente de determinación, Ecuación (13)), a 
(pendiente), y b (ordenada del origen), que relacionan los datos medidos 
localmente con los obtenidos por el modelo aquí propuesto tanto para el 


producto final como para analizar las variables de entrada del método: 


R-RMSE = VRSD? + Me? (12) 
ette E J[o- 0 )r 
r?r= ALE (13) 


donde: 
Ei = datos estimados 
Oi = datos observados 


N = cantidad de datos tenidos en cuenta 
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Resultados y discusión 


Análisis de las variables de entrada 


Con el fin de corroborar el comportamiento de las variables de entrada 
con mayor peso en el cálculo de ETp, se comparó la Rs!, RI!, Ta y Ts 
provenientes de NOAA/CFSV2 y MODIS (GEE) con la base de datos NASA 
Power en toda la RPA, tomando valores de los sitios de las siete estaciones 


mencionadas en la Figura 2. 


En este sentido, se hizo un análisis estadístico (Tabla 4, Figura 3A, 
B, C y D), donde se puede observar que para la variable Rs! (Tabla 4, 
Figura 3A), el análisis estadístico arroja valores de R-RMSE de 36 W m-? 
(del orden de 1.3 mm d! y de alrededor del 10%), y valores positivos de 
mediana, lo que indica una sobreestimación de los datos de NOAA/CFSV2; 
este error sistemático puede ser corregido con un ajuste lineal. El error 
aleatorio (RSD) es de 15.9 W m-*? (alrededor de 0.5 mm d*1); la pendiente 
es cercana a 1 mm d”!, y el R? es de 0.95, marcando una buena relación 


entre las variables. 


Tabla 4. Estadísticos que relacionan los valores de Rs!, RI' Ta y Ts de 
reanálisis (NASA Power) con los datos de satélite (NOAA/CFSV2 y 
MODIS-GEE) para el periodo 2015-2021. 


NOAA/CFSV2/FOR6H MOD1iA2 
Rs'(W m7?) RI'(W m2) Ta (K) Ts (K) 
Mediana 32.8 9,2 0.5 0.2 
RSD 15.9 5.6 0.9 1.3 
l R=RMSE | 36.5 | 108 | 10 | 1.3 | 
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Figura 3. Relación de valores de NASA Power con datos NOAA (A) Rs); 
B) RI'; C) Ta), y MODIS (D) Ts (obtenidos de GEE) para el periodo 
2015-2021. 


En relación con la Rl' (Tabla 4, Figura 3B), el R-RMSE es de 10.8 W 


m-”? (alrededor del 4% - 0.3 mm d?); la mediana negativa indica una leve 
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subestimación de los datos; el RSD equivale a 5.6 W m? (0.1 mm d”?); la 


pendiente es cercana a 1 mm d”!, y el R? es de 0.93. 


Por otro lado, la Ta (Tabla 4, Figura 3C) muestra valores de R-RMSE 
de 1 (K), una mediana positiva, determinando una sobreestimación de los 
datos de satélite, aunque en bajas proporciones. El error aleatorio es de 
0.9 (K); la pendiente es alrededor de 1; el R? es de 0.96, y la distribución 
de los datos se encuentra cercana a la línea 1:1, mostrando una buena 
relación entre los datos de NASA Power, en comparación con los obtenidos 


desde satélite. 


En cuanto a la Ts (Tabla 4, Figura 3D), los valores de los estadísticos 
indican errores bajos (1.3 K), así como la mediana positiva, que 
determina una sobreestimación de los datos de satélite poco notoria, ya 
que los valores se observan cercanos a la línea 1:1, y el valor de R? es de 
0.94. 


Además, con el fin de corroborar estos resultados, se compararon 
las variables de entrada con datos de campo, medidos en una estación 
experimental de cebada (-37.49; -58.90), para el ciclo de cultivo del año 


2019 (sensores y variables detalladas en Faramiñán et a/., 2022). 


En este análisis se obtuvieron resultados similares a los obtenidos 
con los datos de NASA Power, en la RPA en general, para el conjunto de 
datos de Rs! y RI' (Tabla 5, Figura 4A y B). Los errores indican un leve 
aumento en campo tanto para Ta como para Ts (Tabla 5, Figura 4C y D). 
Asimismo, los valores de la pendiente son cercanos a 1 en todos los casos 
y el valor del R? es superior a 0.9 para las cuatro variables. Este análisis 
resuelve que los datos de entrada utilizados para el cálculo de ETp son 


confiables para su aplicación en el método de PT. 
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Tabla 5. Estadísticos que relacionan los valores de Rs!, RI! Ta y Ts 


medidos en campo con los datos de satélite. 


NOAA/CFSV2/FOR6H MOD11A2 
Rs'(W m-?) RI'(W m2) Ta (K) Ts (K) 
Mediana 32.4 -8.8 2.0 1.2 
RSD 15.7 4,3 1.0 1.8 
R-RMSE 36.0 9.8 2.2 2.2 
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Figura 4. Relación de valores de campo con satélite A) Rs:, B) RÍI:, C) 
Ta, D) Ts. 

Si bien se observan algunas diferencias en el análisis de los datos 
de satélite, se puede determinar que muestran una gran potencialidad 
para la utilización de los mismos en el cálculo de la ETp y ETr debido a su 
concordancia con los datos de referencia. 
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Análisis del producto de ETp 


Utilizando la base de datos de la ORA como datos reales (observados), se 
realizó la valoración del producto de ETp obtenido mediante la ecuación 
de PT para las siete estaciones de la RPA (ver Tabla 3 y Figura 2). Para 
esta valoración, se compararon los datos medidos in situ con los obtenidos 
con el modelo a partir de los datos de satélite (Figura 5). Por otro lado, 
en la Tabla 6, se listan los estadísticos obtenidos para esta comparación. 
Los resultados muestran una distribución de la nube de puntos 
consistente, con un R? cercano a 0.9 para toda la región en su conjunto 
como en particular en cada estación (varía 0.7 y 0.9). Por otro lado, la 
pendiente general de esta relación es cercana a 1 y la ordenada del origen 
es de 0.05 mm d*?, 
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y = 1.0915x+0.053 


R? = 0.8601 a A 
N = 2254 


1 2 3 4 5 6 7 8 
ETp (campo) [mm d**] 


o Concordia “UGeneralPico Junín 4Laboulaye XPigué =Rosario + Tandil 


Figura 5. Comparación entre datos de ETp de campo y satélite. Línea 


negra: relación 1:1. Línea punteada: línea de tendencia. 
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Tabla 6. Estadísticos que relacionan la ETp modelada con los datos 


puntuales. 
ETp RPA Concordia | General Pico Junín Laboulaye Pigúé | Rosario | Tandil 
Mediana 
0.3 0.6 0.2 0.0 0.3 0.0 0.1 0.5 
(mm d?) 
RSD 
0.5 0.3 0.8 0.6 0.4 0.6 0.4 0.2 
(mm d*1) 
R-RMSE 
0.5 0.7 0.8 0.6 0.5 0.6 0.4 0.5 
(mm d”!) 
a 1.09 1.17 1.12 1,12 0.98 0.95 1.1 1.17 
b 
0.05 0.23 -0.07 -0.21 0.42 0.19 -0.16 0.11 
(mm d**) 
R? 0.86 0.89 0.79 0.83 0.72 0.87 0.90 0.93 


Analizando a detalle lo referente a la mediana (0.3 mm d!), este 


factor otorga un valor positivo para toda la región, mostrando una leve 


sobreestimación de los datos, determinando un error sistemático casi 


imperceptible, como lo es en la mayoría de las estaciones, en particular 


en Junín, Pigué y Rosario; es mayor en las zonas de Tandil y Concordia. 


Con respecto al error aleatorio, determinado por el RSD, se obtienen 


valores de 0.5 mm d'* para la RPA en general, siendo mayor en las 


estaciones ubicadas en la zona semiárida (Pigúé y General Pico) y menor 


en la subhúmeda/húmeda. En cuanto al R-RMSE, su valor es de 0.5 mm 


di en la RPA en general, que se encuentra por debajo del error de 


estimación del método; según la Tabla 1, para estimaciones realizadas 


con satélite, el error típico es entre el 10 y 20% (Allen et a/., 2011). 
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Los errores asociados con este producto de ETp son comparables 
con los obtenidos por otros autores que calculan ET con el método de PT 
y PM (Monteith € Unsworth, 1990) en distintas regiones del planeta. En 
particular, Laipelt et a/. (2021), que obtuvieron la ET con el método SEBAL 
(Surface Energy Balance Algorithm for Land), utilizando datos LANDSAT 
y ERA5 con GEE, arrojando errores medios de 0.67 mm día”? y valores R? 
superiores a 0.4 para la mayoría de los sitios evaluados en distintas zonas 


de Brasil. 


Por su parte, Yao et al. (2017) aplicaron la ecuación de PT para 
calcular ETp en varios biomas y climas del noreste de China. Utilizaron 
como datos de entrada distintos productos MODIS, logrando una 
resolución espacial de 1 km y temporal de ocho días. Estos autores 
observaron errores menores en zonas con cobertura vegetal continua y 
con disponibilidad de agua (alrededor de 0.6 mm d**). Mientras que en 
regiones sin cobertura o con cobertura discontinua y con clima 
semiárido/árida, el algoritmo muestra errores mayores (alrededor de 
0.8/0.9 mm d”). 


Por otro lado, Degano et al. (2021c) evaluaron el producto 
MOD16A2 (ETp) para la RPA (Mu, Heinsch, Zhao, 8 Running, 2007; 
Running, Mu, 8 Zhao, 2017). Este producto, disponible globalmente, cuya 
resolución espacial es de 1 km (discontinuado) y 500 m, y temporal de 
ocho días, está basado en el modelo de PM. Los autores encontraron que 
para esta zona de estudio, el producto MOD16A2 presenta un error 
sistemático, indicando las debilidades del mismo para ser utilizado de 


forma directa en la RPA. Por lo tanto, desarrollaron una ecuación de ajuste 
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lineal para corregir ese producto, logrando disminuir los errores asociados 
(de 3.4 mm d* a 0.6 mm d!). 


En la Figura 6 se muestran los mapas anuales de ETp para la RPA 
(2016-2021). La ET aumenta con sentido sudeste-noroeste, con un patrón 
similar entre años. En la zona semiárida, donde hay mayor demanda de 
agua de la atmósfera, la ETp es mayor; en las zonas más húmedas, los 


valores de ETp disminuyen. 
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Figura 6. Variabilidad espacial de la ETp en la RPA. 
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En este estudio se tuvieron en cuenta seis años, de los cuales se 
encuentran años extremos y promedios en cuanto a la precipitación (800- 
1 100 mm) y humedad del suelo (40-70%). En particular, el 2016 y 2017 
fueron los más húmedos, el 2018 y 2019 con precipitaciones promedio, y 
los años restantes más secos. El análisis temporal muestra un aumento 
de la ETp desde el 2016 (mayores precipitaciones) al 2021 (menores 
precipitaciones), es decir, a mayor demanda de agua de la atmósfera se 
observan valores de ETp más altos, incrementándose el límite de los 
valores altos de este parámetro, mostrando un corrimiento del mismo 


hacia el sudeste para el año 2021. 


Este análisis multitemporal de la ETp puede ser determinante para 
los productores agropecuarios, dado que la disponibilidad de agua del 
sistema es vital para el desarrollo de los cultivos. En dicho sentido, y 
teniendo en cuenta esta variable, las personas involucradas en la 
producción de los diferentes cultivos serán posibles usuarios de este 
producto de ETp a escala regional (en el marco del proyecto JOVIN se 
está realizando la transferencia de la aplicación de este producto a la ORA, 
perteneciente al Ministerio de Agricultura, Ganadería y Pesca-MAGyP (03- 
JOVIN-711)). 


Análisis del producto de ETr 


En la comparación entre datos de satélite con datos in situ respecto de la 
ETr (Figura 7, Tabla 7), se puede observar una mayor dispersión en la 


nube de puntos, indicado por el valor de R? de 0.54, disminuyendo en las 
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estaciones de Concordia y Tandil; en el resto de las estaciones varía entre 
0.6 y 0.7. 


y = 0.9463x - 0.0629 
R* = 0.5418 
N = 2143 


ETr (campo) [mm d*] 


> Concordia mGeneral Pico <Junin a Laboulaye xPigué -—Rosario +Tandil 


Figura 7. Comparación entre datos de ETr calculados con balance 
hidrológico (ORA) y los de satélite. Línea negra: relación 1:1. Línea 


punteada: línea de tendencia. 
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Tabla 7. Estadísticos que relacionan la ETr modelada con los datos 


puntuales. 
ETr RPA | Concordia | General Pico | Junín | Laboulaye | Pigiié | Rosario | Tandil 
Mediana 
-0.2 0.6 -0.5 -0.4 -0.4 -0.2 -0.3 0.0 
(mm d”*) 
RSD 
0.5 0.5 0.4 0.4 0.4 0.3 0.5 0.4 
(mm d*1) 
R-RMSE 
0.6 0.8 0.6 0.6 0.5 0.4 0.6 0.4 
(mm d”!) 
a 0.95 1.31 0.87 0.97 0.81 0.86 1.18 0.73 
b 
-0.06 0.05 -0..21 -0.36 -0.01 0.03 -0.53 0.29 
(mm d**) 
R2 0.54 0.54 0.6 0.69 0.64 0.63 0.61 0.44 


El error sistemático, representado por la mediana, muestra un valor 
negativo (para el total de los datos de la RPA), determinando una leve 
subestimación por parte del producto de satélite. Esto sucede en la 
mayoría de las estaciones, excepto en Tandil, que es cero, mostrando la 
nube de datos sobre la línea 1:1, y en Concordia, que obtiene un valor 
positivo. 


El RSD, al igual que en la ETp, toma un valor de 0.5 mm d*! en 
general; en particular, este error aleatorio varía entre 0.3 y 0.5 mm d*!, 
Por otro lado, el R-RMSE alcanza valores de 0.6 mm d*? para toda la RPA, 


siendo mayor en zonas con climas más húmedos (Concordia). 


En el contexto de la ETr, diferentes autores obtuvieron este 


parámetro integrando las propiedades del suelo, así como la Hs. En su 
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trabajo, Walker y Venturini (2019b) calcularon ETr, siguiendo la 
metodología de Walker, García, Venturini y Carrasco (2019a), tomando 
como datos de entrada variables de reanálisis disponibles en GEE (Global 
Land Data Assimilation System version 2.1-GLDAS y National Centers for 
Environmental Prediction Climate Forecast System version 2-NCEP). 
Analizaron su método tomando distintas estaciones de la red FLUXNET 
alrededor del mundo. La valoración del método otorgó errores del orden 


de los 3 mm d*!, 


Por su parte, Degano et al. (2021c) evaluaron el producto MOD16A2 
en la RPA para distintos ciclos de cultivo de soja (entre 2009-2010 a 2017- 
2018). Para la ETr determinaron que los valores de los errores resultan 
entre 0.7 mm di y 2 mm d*! en la zona subhúmeda/húmeda, mientras 
que en la zona semiárida los valores disminuyen a 0.5 mm d*!. El 


promedio general del error en la zona es de 1.1 mm d”?, 


Este análisis indica que el producto desarrollado en este trabajo 
muestra errores competentes en relación con los obtenidos por otros 


autores, variando entre 0.4 y 0.8 mm d”!, 


La Figura 8 muestra la variabilidad de la ETr en la RPA en los 
diferentes años (2016-2021), donde la variación espacial está controlada 
por la disponibilidad de humedad en el perfil de suelo. El comportamiento 
espacio-temporal de la ETr está influenciada por la distribución de las 
precipitaciones anuales con un fuerte control del tipo de suelo 


(dependiente de sus características texturales y estructurales). 
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En particular, cuando se trata de años más húmedos, la tendencia 
equivale a una menor variación espacial de la ETr (en general, valores 
más altos para el conjunto de estaciones). En cuanto a los años con déficit 
hídrico, en general se presenta una mayor variación espacial y con mayor 


impacto en zonas de suelos con menor capacidad de retención de agua. 


Analizando el cambio multitemporal de la ETr en la zona del río 
Paraná (noreste), se puede observar con claridad una disminución del 
proceso de pérdida de agua del sistema a través de los años. Esta 
variación puede atribuirse a las sequías observadas por los efectos sobre 


las precipitaciones de los eventos de La Niña sobre la región. 


En concordancia con la ETp, y haciendo uso del producto de ETr, 
esta herramienta puede ser de gran ayuda para los productores agrícolas 
a la hora de planificar la siembra, así como la distribución de los cultivos 
en la RPA. Teniendo en cuenta la disponibilidad de agua tanto a nivel 
atmósfera como a nivel del suelo, y dependiendo de las predicciones a 
futuro (año Niña-seco o Niño-húmedo), puede influir en la toma de 


decisiones para el sector agropecuario. 
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Conclusiones 


En este trabajo "se propone un método para determinar 
evapotranspiración potencial y evapotranspiración real a nivel regional 
con datos de satélite y de reanálisis con la herramienta Google Earth 
Engine (GEE), con una resolución temporal de ocho días y espacial de 250 
x 250 m, utilizando el método de Priestley-Taylor. Este producto fue 
valorado con datos de campo, provistos por la Oficina de Riesgo 


Agropecuario, en la región pampeana argentina. 


En primera instancia, se evaluaron los datos de entrada de satélite 
con datos NASA Power y con valores medidos localmente. En este sentido, 
se concluye que tanto los datos de NOAA/CFSV2 como los de MODIS 
(ambos obtenidos de GEE) muestran una gran potencialidad para realizar 
el cálculo de ETp con el método de PT debido a su concordancia con los 


datos de referencia. 


Respecto de la evapotranspiración potencial, en general, el producto 
arrojó errores de 0.5 mm d*, mostrando una leve sobreestimación, 
determinado por la mediana. En el análisis multitemporal realizado se 
observa una distribución de esta variable, que aumenta en sentido 


sudeste-noroeste en toda la región pampeana argentina. 


El análisis estadístico relacionado con la evapotranspiración real 
mostró errores del orden de 0.6 mm d“! y una breve subestimación del 
producto de satélite. Por su parte, el análisis temporal de este parámetro 


en la región muestra una variabilidad, aumentando de sudoeste a noreste. 


En general, debido a que los errores que presenta el producto están 


dentro del error asociado con el método, se puede concluir que este 
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producto puede utilizarse de manera directa para obtener los valores de 
evapotranspiración potencial y real en una zona de llanuras con 
características subhúmedas/húmedas a semiáridas, y que puede 
emplearse para analizar la variabilidad espacio-temporal de la ETp y ETr 
en diferentes condiciones extremas (año Niño/ año Niña), así como 


realizar estudios ambientales aplicados. 


Además, el uso de la herramienta Google Earth Engine logra un 
avance importante en el procesamiento y análisis de las variables 
ambientales/hidrológicas de satélite y reanálisis, pues tiene una amplia 
variedad de catálogos de datos disponibles, con buena precisión e 
integración de los datos para lograr un producto aplicable, con técnicas 
sencillas de programación, lo cual agiliza el procesamiento de los datos 


de entrada y la obtención de resultados. 
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